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5^ ■ Abstract 



We show that the resummation of large perturbative corrections in QCD leads to am- 
biguities in high energy cross sections that are suppressed by powers of large momentum 
5^ ' scales. These ambiguities are caused by infrared renormalons, which are a general feature 



of resummed hard-scattering functions in perturbative QCD, even though these functions 
are infrared safe order-by-order in perturbation theory. As in the case of the operator prod- 
uct expansion, the contributions of infrared renormalons to coefficient functions may be 
absorbed into the definition of higher-dimensional operators, which induce nonperturbative 
corrections that are power-suppressed at high energies. The strength of the suppression is 
determined by the location of the dominant infrared renormalon, which may be identified 
explicitly in the resummed series. In contrast to the operator product expansion, however, 
the relevant operators in factorized hadron-hadron scattering and jet cross sections are gen- 
erally nonlocal in QCD, although they may be expressed as local operators in an effective 
theory for eikonalized quarks. In this context, we verify and interpret the presence of 1/Q 
corrections to the inclusive Drell-Yan cross section with Q the pair mass. In a similar 
manner, we find exp(— 6^ InQ) corrections in the impact parameter space of the transverse 
momentum distributions of the Drell-Yan process and e'^e" annihilation. We also show 
that the dominant nonperturbative corrections to cone-based jet cross sections behave as 
1/{Q5), with 6 the opening angle of the jet and Q the center of mass energy. 



1. Introduction 

It is almost axiomatic that perturbation theory alone cannot give a full description of QCD. 
Luckily, however, perturbative calculations may be supplemented by nonperturbative functions 
and parameters to derive physical predictions. While other sources of nonperturbative behavior 
are possible, it is attractive to identify those required to make perturbation theory well defined. 
In this manner, we may use perturbation theory as a diagonstic tool to identify a minimal set of 
nonperturbative parameters. This is the viewpoint that we shall advocate and exploit below. It is 
of particular interest for infrared safe quantities, which are finite order-by-order in perturbation 
theory. For such quantities, we expect perturbation theory to predict leading behavior as an 
asymptotic expansion in the coupling, and for the failure of perturbation theory at high orders 
to be particularly enlightening. 

To be specific, we shall work in the class of quantities for which leading logarithmic behav- 
ior in energy is described by Sudakov resummation |l|, |^, |^, ^, |^. This resummation may be 
formulated in terms of path-ordered exponentials, or Wilson lines. Path-ordered exponentials 
have been extensively studied in perturbation theory ||^, ^ and are also of special interest for a 
nonperturbative description of the theory. We will see below that they naturally relate so-called 
infrared renormalons [^ ||, |^ to nonlocal operator expectation values [|1^] for cross sections of 
lepton pair production ||lT], |l^ and energy-energy correlations @ in e"'"e~ annihilation. Gener- 



alizations to other quantities, including jet cross sections, and to nonleading logarithms are also 
discussed. 



Already in [13| it was shown that resummation of large perturbative corrections to the in- 
clusive lepton pair cross section leads to an ambiguity in perturbation theory at the level of 
A/Q, with Q the pair mass and A the QCD scale, compared to leading power. This paper is in 
large part an attempt to further understand such power-suppressed effects^. We shall identify 
the corrections implied by these ambiguities with specific nonlocal operators below. In a sim- 
ilar manner, we shall find exp(— fo^lnQ) corrections (previously-discussed in 0) in the impact 
parameter space for the transverse momentum distribution in the Drell-Yan cross section and in 
e"'"e~ annihilation. In addition, we shall see that nonperturbative corrections to cone-based jet 
cross sections may be expected to begin at 1/{Q5), with 5 the opening angle of the jet and Q 
the center of mass energy. An application to inclusive B meson decay has been given in ||15|| . 

Unaided, perturbation theory fails in a number of ways. At low orders, infrared divergences 
and long-distance dependence at fixed orders in perturbation theory can be factorized into non- 



perturbative functions like parton distributions, which are determined by experiment ||T6[. In 
addition, QCD perturbation theory is (probably) at best an asymptotic expansion for Green 
functions, even off the mass shell. Relevant for us are individual diagrams that behave at n loops 
as a"n!. A series that behaves in this fashion is not Borel summable, let alone convergent, be- 
cause such terms produce singularities in the Borel transform that prevent it from being inverted 
in a unique fashion. Singularities that are associated with individual diagrams are often termed 
infrared or ultraviolet renormalons. In QCD, it is the infrared (IR) renormalons that directly 
interfere with the inversion of the Borel transform 0, ^ . 

A lack of Borel summability does not imply that a series is meaningless, however. Rather, it 
may represent an asymptotic expansion for a set of functions that differ only in ways that do not 
affect this expansion. For instance, they may differ in powers of exp[— 1/a], with a the expansion 
parameter. Different members of this class of functions correspond to various definitions of the 

^For recent related studies see Ref. fH. 



inverse Borel transform. It is this sort of ambiguity with which we shall deal below, showing 
how its resolution requires the introduction of new, nonperturbative, parameters associated with 
vacuum expectations for nonlocal operators^. 

The paper is organized as follows. In Sect. 2 we review how infrared renormalons arise in 
resummed cross sections. Analyzing their contributions to a Wilson line expectation value, we 
then find a general form for power suppressed corrections associated with infrared renormalons in 
the transverse momentum distribution of lepton pairs. In Sect. 3 we relate the contribution of the 
leading infrared renormalon to the matrix element of a nonlocal gluon field operator, which can 
also be represented as a local operator in an effective field theory for eikonal quarks. In Sects. 4 
and 5 we consider power corrections to jet cross sections and inclusive lepton pair production, 
and identify in each case a leading infrared renormalon that results in A/Q corrections. Sect. 6 
contains concluding remarks. 

2. Infrared renormalons in a resummed cross section 

In many applications of perturbative QCD to high energy scattering, it is desirable to sum finite 
corrections that remain after the cancellation of infrared divergences. Examples are found in 



energy-energy correlations in e^e~ annihilation [^ |T8|, the transverse momentum distribution of 
lepton pairs |jl9|, |Tl|, Sudakov effects in elastic scattering ||2^, the total cross section for lepton 
pair production |21, 22, 23, 2|, 13| and for heavy quark production |^|, and also in heavy quark 



effective theory [^, ITsf. As we shall see, such resummations of perturbative corrections generally 
imply the presence of nonperturbative corrections suppressed by powers of kinematic invariants. 
Let us first consider the transverse momentum distribution of lepton pairs in the Drell-Yan 
process, h + h' ^ i~^i~ + X. The relevant formalism is virtually identical for energy-energy corre- 
lations in e"'"e~ annihilation, closely related to jet structure 0. For large transverse momentum 
of the lepton pair, Qt ~ Q, the differential cross-section for the process may be written to leading 
power in the lepton invariant mass Q as 

where ctq = Air'^a'^/dQ'^s is the inclusive Born cross section for a quark of unit charge, where s 
is the squared invariant mass of the incoming hadrons, and where the function E depends on 
parton distributions fa/h{.x,jj?) through the factorization formula 



nQ\Ql) = t— — fa/H{Xa,Q')hMx,,Q')Hj-^,%) . (2) 

JO Xa Xb \XaXbS Q^ J 

The bulk of the lepton-pair cross section, however, is in the region Q\ <^ Q^ . In this region 
the function S contains "double-logarithmic" corrections like l/Q^ x a" ln^"~^(Q^/(5^). The 
resummation of these large corrections is most easily expressed in the Fourier b-space conjugate 
to Qi, 

S(Q^Q?) = /^e^^*■^S(g^6^), (3) 



^An application of similar ideas to QED is given in [[l7[ . 



where 6^ = b^. Then, the double-logarithmic corrections to S((5^, &^) are of the form a" In " Q^fc^ 
and exponentiate. At leading power in Q^, the expression for E((5^, 6^) is given by 0, |TlI 



)2 i?\ = r ^^ f „ (^ ri}?\ f. „.(rr., nll?\^(^ _ n2 u ^, ^.^ ^-5(0^6^) 



S(Q^fe^)= / — — /aM(x.,C/6^)^/,,(a;,,C/6^)5(l-Q7x.x,.)e-^^^'''^^ (4) 

Jo Xq Xb ' 



where e ^'^^ '^ Ms a resummation factor that includes the large logarithms in Q'^h 
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In --2-rcusp(as(/i;t^)) + T{as{kl)) 



C = Ae-'^ , (5) 



with 7 the Euler constant. Here, Fcusplas) is the universal "cusp" anomalous dimension |^, |[, 
while T{as), which contributes to nonleading behavior in InQ, depends on the specific process 
and choice of the constant C. With the choice above, both anomalous dimensions begin at order 

rcusp(«.) = -Cf + 0{al) , r{as) = -^-Cf + 0{a^) , (6) 

with Cp = 4/3. A short calculation, using the one-loop form of the running coupling, 

"^^^'^ = AM^fM' /5i = (ll-2n;/3)/(4vr), (7) 

shows that the leading (double) logarithms of Q'^b'^ found order-by-order are replaced in the 
resummed form by 



SiQ^b^'- ^^ 



^^^^'/^')^^]^^-^^^(^'''/^) 



where we have used the one-loop expression (^ for r^usp and have neglected nonleading loga- 
rithmic terms due to F in (|^). At large b, the Fourier transform S(Q^, 6^) is strongly suppressed 
by the resummed exponent, which makes the structure function J]{Q'^,Q^) amenable to pertur- 
bative treatment. We also note, however, that the exponent diverges at 6^ = C/A^, because of 
the singularity in the running coupling (|^) at kf = h?. It is the consequences of this singularity 
for power corrections in b^ that we shall discuss below. First of all, we notice that approaching 
the "critical" value of the impact parameter b ~ 1/A, one has to take into account power correc- 
tions (9(A^6^), which are negligible with respect to logarithmic corrections ln(A^6^) for small 6^, 
but which become important for large 6^. We stress that the factorization in (^) was found in 
the leading 1/Q^ limit and therefore in principle includes power corrections of 0{K%'^) provided 
that b'^Q^ ^ 1, or equivalently Ql -C Q^. At the same time, the expression (|]) was obtained 
by summing all logarithmic corrections and neglecting power corrections. Thus, to study the 
6^ power corrections we have to find an approximation for the exponent S{Q^, b^) that includes 
both logarithmic and power corrections in 6^. To this end, we recall the origin of e"'^'''^ '^ ^ in the 
factorized expression (^). 

In the Drell-Yan process, the resummation factor, e"'^'^'^ '^ \ appears as the contribution of 
real and virtual soft gluons interacting with the quark and antiquark that annihilate to produce 
a virtual photon. The interaction of soft gluons with quarks can be treated using the eikonal 
approximation and can be summarized in terms of path-ordered exponentials (Wilson lines) ||2^ 



as follows. Let us introduce a notation for "straight" path-ordered exponentials in the direction 
of momentum p, 

rs 

'So 



^p[Sp + X, SqP + x] = P exp i — ig dsp^A^(ps + x) j , (9) 



where P denotes ordering in group indices of the gluon field A^{x). Similarly, for products of two 
such ordered exponentials, defined by the light-like momenta pi and p2 of the incoming quark 
and antiquark, respectively, we define 



Uy)y{x) 



$T (x, -OO) $pj(x, -OO) 



^^p^{-oo,x) $p,(x, -oo) 



pI=pI = o^ 



(10) 

where T is time-ordering. Finally, in terms of such path-ordered exponentials, we define the 
eikonalised Drell-Yan cross-section as 



iyDY(&,/i) = (0| UlyiO) Uj^yib) |0) , b-' = (0,0, b) (11) 

with bub'^ = —b'^. If we insert a complete set of states, 1 = J2x l-^)(^l between the operators U 
and U\ 

WnYib,fi) = Y. 1(01 UnriO) |X)|' exp(-2k, ■ b) , (12) 

X 

the Wilson line expectation value PVdy describes the cross section for the emission of gluons of any 
energy by two eikonahzed quarks, weighted by a factor exp(— -ikt -b), with kj the total transverse 
momentum of gluons in the final state \X). Such a cross section is ultraviolet divergent, and 
requires renormalization, which makes Wdy a function of renormalization scale /i. As discussed 
in [^, for instance, Wdy takes into account the contribution of soft gluons to the Drell-Yan 
cross section. The value of /x is arbitrary but the most convenient choice is /i ~ Q. Let us show 
that for /i = Q the expression (0) indeed coincides with the resummation factor (H) up to power 
corrections in 6^. 

Expanding (pH]) in powers of the gauge field, we get the following unrenormalized expression 



for Wdy to the lowest order of perturbation theory: 



<^ = 1 + ,V-^C^ f:^2n6^{e)(^-^]\l-e-^^^-^), (13) 



where /i is the parameter of dimensional regularization, with D = A — 26, d^k = dk^dk^d^^^kt, 
and ki is the {D — 2)— dimensional transverse momentum of the gluon. For massless quarks, the 
/c* integral contains a collinear divergence. This divergence is eliminated, however, when one 
evaluates the logarithmic derivative dlnM/Dv/dlnQ^. Using the one-loop expression W^y and 
the exponentiation of infrared divergences in hard-scattering processes [^, ^, ^ we get 

dlnl^DY _ 4_g f d^-^fcf «w ,i,^.b ^ . . 

To take into account higher order corrections to (|1^) we choose the argument of the coupling 
constant to be the transverse momentum, a^ = as{k^), and identify CfQ-s/'^ as the lowest term 
in the expansion of the cusp anomalous dimension, rcusp(as(^?)) [i, H, 



,n.Q^ '^^1 (2.)^-^ -'-•'■ [e-^-^ - I] . (15) 



Notice that after integration over large k'f this expression contains uhraviolet poles in e which 
need to be subtracted. The resulting renormalized expression contains both logarithmic and 
power corrections in 6^ to the eikonalized cross section, due to soft gluon emissions from the 
incoming quarks. Neglecting power corrections to (piSf ), we find by comparison with eq. (^), 
that the behavior of In VFdy in l^Q in (^) is the same as the behavior of S{Q, b) in InQ up to 
nonleading logarithmic terms due to the "collinear" anomalous dimension F in (H). Identification 
of the resummation factor exp{—S{Q^, 6^)) as a Wilson line expectation value ( jTSD allows us now 
to analyze the perturbative contribution of soft gluons with small transverse momentum kf to 
the Drell-Yan cross section. 

In the evolution equation ([T5|) and its solution VFdy "we find divergences at k^ = h?. Notice 
that ( P!3[ ) was found after resummation of soft gluons to all orders of perturbation theory, and that 
this singularity can easily be translated into the high-order behavior of the original perturbative 
series before resummation. We can reconstruct the expansion in as{Q'^) by using the one (or 
higher) loop running coupling, as(fc^) = as(Q^)/[l — /5ias(Q^) InQ^/fc^]. Such an expansion 
generates a sum of integrals of the form /q da;ln"(l/x) = r(n-|-l), and we find that the singularity 
manifests itself in the perturbative expansion in as{Q'^) as an infrared renormalon in the exponent 



T3| , behaving at nth loop order as a;"/?"n!. 

Let us now study in more detail how IR renormalons appear in the evolution equation (p^), 
including renormalization. After substitution of the relation as{k^) = J^ da{k^ / A"^) " ^ into 



(jTj), and after integration over transverse momenta, we get 



dlnWV _ ^CV(^^.,.). /°° f" r(l-.a-.)/A^y^. ^ j^gj 



dlng2 vr ' "^ Wo a(3i+e T{l + a(3i) V 4 

We notice that {K%'^ /A)"^^ = ^-'^1°-^ with a^ = as(4/6^) and identify the right-hand side of the 
last relation as the Borel representation, 

/•OO 

7r(a,) = / d(Tn{(r) e^'/^^ , (17) 

Jo 

oi n(as) = din Wov/dlnQ^, with 

vr opx + e 1 (1 + (Tpij 

Let us analyze the singularities of the Borel transform 7r(cr) for real positive a. In the limit £ -^ 
the function 7r(cr) has a pole for a = 0. Since small values of the Borel parameter correspond to 
small coupling constants asikW or equivalently to large transverse momenta fc^, the singularity 
at 0" = is ultraviolet. Away from a = we put e = and find that the function 7f(o") has 
singularities generated by F— functions at a* = l//3i, 2//5i, 3//3i, . . . , the infrared renormalons. 
To define the integral over large a we have to fix the prescription for integration of the IR 
renormalons in ([TBI) . Different prescriptions lead to results that differ from each other in powers 
of 6^A^, beginning at 0{\?t^\ the contribution of the "leading" IR renormalon, at a* = l//?i. 

To get a "perturbative" approximation to VTdy; we treat cr as a small parameter and then 
integrate (|TBp over < o" < 1/2/?!. For small a we replace the ratio of F— functions in (|TB|) by 
g27cr/3i+7e^ use the expansion l/(e + o"/?i) = !/£: — ajSi/e^ + . . . and change variables to x = cr/3i 
to get 

dlnVTpT Cf ^, ,„ (27r/iV/C)^ /-i/a 



dlnQ2 vr/3i 



U-r ^\j^ I dxx-\K%'/Cr^ (19) 

n=l ^ -^0 



with C given in eq. (§). Subtracting poles in the MS scheme, we find after summation the 
following result, for /x^ = 2Q^ and 1/6^ > A^, 



dlnQ2 7r/?i yi/in(Q2/A^) a; vr/^i \n{Q'^ / JS?) 

One easily checks that the resummation factor (||) satisfies this relation up to nonleading loga- 
rithms in Q, so that to this approximation, VFpt = exp (— ^(Q^, &^)). 

There is no unique way of defining the kt integral in eq. (|15D or the a integral in eq. ([I6|) . As 
discussed above, however, we may consider the singularities of the running coupling as inducing an 
ambiguity in the integral, which is to be eliminated by adding new, nonperturbative parameters 
to the theory [R, EBf. This process is not utterly arbitrary, because the singularity appears only 



at small values of fc^, where, unless h is very large (ie., of order 1/A), the integrand is suppressed. 
This means that we may consistently define the full integral in the right-hand side of (|TB|) as a 
power series expansion in 6, starting with a "perturbative" contribution S'pt(&) = — lnP^px(&), 
eq. (ID, at order {h^f: 

^-siQ^b-) ^ ^^^(^) ^ g^p (-5pt(&) - b'S^iQ) - b'S,{Q) + ...), (21) 

where the equality is approximate, since we have neglected nonleading collinear corrections to 
S{Q'^, b"^). The fact that the Borel transform ([18|) has IR renormalon singularities at a* = n/j3i for 



n = 1,2,... implies that nonperturbative effects should contribute (9(6^" A^") power corrections 
to the exponent in (0), or equivalently to the functions S2n{Q)- Moreover, the explicit form of 
the evolution equation (|T5|) implies that in the large Q limit these functions behave as 

^2n(Q)~^2„lnQ + 52„. (22) 



This is because the IR renormalons arise from small values of kt in (jT^), where we may take 
e = 0. Then the contribution of the IR renormalons in (^) to dH^Dv/dQ^ is independent of Q, 
which, along with (pT]), implies the functional dependence (p^). 

3. Operator content of the leading IR renormalon 

We now identify the operator content of the leading infrared renormalon, a* = l//3i, which 
contributes to the 6^A^ correction to the right-hand side of (|15]), by recognizing it as the first 
term in the expansion of IuM/dy with respect to b^. We denote this contribution as 



S2{Q) = -^lnWj,Y{b) 
Applying this derivative to the logarithm of Wdy, eq. (0), evaluated at 6 = 0, we find 



(23) 



S2{Q) = \{0\ 4y(0) m' UMO) |0) , (24) 

where the derivative d = (9i,92) acts in the transverse 6— space. Finally, using the definition 
([10|) of Uby we get a form reminiscent of (plSf ), 



'^2(g) = ^(o| 



$t^(0, -oo)(/,,(0) - :fI^{0))%M -oo) 



|0) , (25) 



where J^p{x) = {iD°'^p{x, — cxd))$^(x, — oo) is a nonlocal "eikonalized" field strength, which may 
also be written in terms of the local field strength F'^'^ and the ordered exponential $p, eq. (§), 
as 

/■O 

J^p(x) = -ig ds ^p(x,x + sp) p^F''°'{sp + x) ^-p(x + ps,x) . (26) 

The perturbative contributions to 5*2 are scaleless integrals and vanish in dimensional regular- 
ization. This is similar to the perturbative renormalization of the divergent matrix elements of 
a local operator such as (0|F^|0). Similarly, we may take VFdy(O) as unity. 

From (^Sf ) and (p6l), we see that the leading nonperturbative correction is described by the 
matrix element of a nonlocal operator of the gauge field. This is in contrast with the total e"'"e~ 
cross section, in which the leading IR renormalon appears at a = 2//?i and gives rise to power 
corrections described by the local operator (0|F^|0). Nevertheless, it is possible to represent 
(|25| ) as a local operator in an effective field theory. In this model, the quark with light-cone 



momentum p is described by the effective field, qp{x), and the interaction with gauge fields, 
Afj_{x), is organized to reproduce the quark-gluon interaction in the eikonal approximation. The 
Lagrangian of the model is given by 

^ = Qpix){p ■ iD)qp{x) , (27) 

where g^ = qj,, with p the quark momentum, and where Z)^ = d^ — igA^{x) is the covariant 
derivative in the quark representation. For a heavy quark with mass M, velocity f ^ and momen- 
tum p = Mv, this is the Lagrangian of heavy quark effective theory. The solution of the equation 
of motion for the g— field reproduces the eikonal phases of the quark fields 

qp{x) = $p(x, -oo) ttp , (28) 

where ap = qp{—oo) is the creation operator of the eikonalized quark. Using this nonabelian 
"Bloch-Nordsieck" model, we may rewrite Udy{x), eq. (p!0|), as a composite local operator in 
the effective theory, Udy{x) ~ qp2{x)qpj^{x), where fields g^^ and qp^ describe a quark and an 
antiquark with momenta pi and p2, respectively. The fields q carry color indices of the quarks, 
and the summation over these indices is implied. In these terms, the matrix element in (^4]) is 
given by 

E2{Q) = \mqr>,qp,) {^^f (g^^g^JIO) , (29) 

provided that (0|aj,^apjO) = {^\a\i^ap^\^) = 1 . As was pointed out above, the value of this 
nonperturbative matrix element will depend on how we define perturbation theory. One definition 
is to apply a principal value prescription for integration around the poles in the inverse Borel 
transform, eq. (|TBp |TB[ . 

A different method is implicit in the prescription given by Collins and Soper [§], in which the 
impact parameter h is replaced by the function 

"•^^TW^- *~2GeV. (30) 

in the expression for the exponent S^Q"^, 6^) and for the quark distribution functions fa/h{x, C/b'^) 
in eq. (^), while an additional nonperturbative factor 

exp(-^l(6^r)-^2(&^)ln^) (31) 



is introduced into the exponent S in eq. (^ to suppress the large-6 region in the Fourier integral in 
eq. (^). In the large-fe hmit, 6* approaches a maximum value of 1/Qo and the perturbative kf inte- 
gral in (|^) becomes well-defined. For moderate 1/Q ^ fe ^ 1/A, the following parameterization 
was also proposed, again consistent with the result of eq. (|2ll) above, 

g^{b^T) = g,b' + 0{b'), g,{b') = g,b' + 0{b'). (32) 



The values of gi and g2 cannot be calculated in perturbative QCD, but have been estimated ||T2[ 
from experiment^, 

gi = 0.15 (GeV)^ ^2 = 0.40 (GeV)^ . 

In summary, the results (^) and (^) of our analysis are in agreement with the Collins-Soper 
parameterization of nonperturbative effects in (|3T|) and (p2D. In addition, our considerations give 
a field-theoretic interpretation to the coefficient g2, as 



with S2 and A2 given in (pSf ) and ([2^). Of course, the specific value of g2 quoted above depends 
on the Collins-Soper method, ( pOD and (|3TD, of defining the perturbative content of resumma- 
tion. Other values will give different g2^s in general, an ambiguity that we expect from general 
considerations. 

It is now useful to discuss the nature of our approximations. We have taken into account the 
contributions of soft gluon radiation from active quarks, and have ignored the contributions of 
collinear gluons as nonleading in InQ^. In perturbation theory, this property follows from the 
expression @ for the resummation factor, in which purely collinear interactions contribute to the 
anomalous dimensions F, while soft emissions contribute to both F and Fcusp- Generalizing this 



observation to the power corrections in the exponent S, eq. (|2l|) , we find that collinear interactions 
modify the coefficients i?2n in (P^), but not the leading coefficients A2n- In addition, considering 
the expression (H) for S((5^, 6^), we notice that for large b^ the normalization scale in the definition 
of the quark distribution functions becomes small, which should produce nonperturbative power 
corrections to fa/h{xa,C/b'^). Such corrections are naturally attributed to "intrinsic" motion of 
quarks inside hadrons. This motion affects the transverse momentum of the lepton pair, and 
may be parameterized phenomeno logically by the function giilP'^r) in (0). This implies that 
the coefficients i?2n in (^D refiect the intrinsic momenta of quarks in hadrons 0. In realistic 
hadronic states, we expect quarks to have a "residual" virtuality, which should replace 1/6^ as a 
regulator of collinear divergences in S{Q'^, 6^), (^), for large b^. At the same time, the coefficients 
A2n are free from collinear divergences, i.e., are infrared safe, and are thus independent of quark 
virtuality. Therefore, we expect the coefficients A2n to be independent of the incoming hadrons. 
This explains why only the vacuum expectation values of operators appear in relations ( ]TT| ) and 

It is clear that this pattern will recur whenever a cross section can be written in terms 
of a resummed expression like eq. (^, that involves the integral over the scale of the running 
coupling. In the following two sections, we shall find two other applications, jet cross sections 
and the inclusive dileption cross section, which show a similar pattern, and which illustrate the 
dependence of the implied nonperturative corrections on the process. 

■^Another parameterization was proposed in P9| . 
Similar phenomena occur in the behavior of the heavy quark distribution function in the end-point region 



4. Power corrections in jet cross sections 

As another application of our infrared renormalon analysis, we examine power corrections to 
infrared-safe jet cross sections. Let us consider the cross section for e"'"e~ -^ 2 jets, defined by 
calorimeters in the form of back-to-back cones of half-angle 6 in the overall center of mass. Let 
^/s = Q he the center of mass energy, and let xQ be the total energy flowing into the two jets. 
For X — * 1, the cross-section dcr2j/dx receives large perturbative corrections associated with the 
emission of soft gluons outside the jet cone, with total energy (1 — x)Q. 

As in the Drell-Yan process, the asymptotic behavior of the cross section dcr2j/da; for x — > 1 
is best seen in terms of its Fourier transform with respect to the total energy (1 — x)Q carried 
by soft gluons, 

da,,iil-x)Q,6) ^r-dyo ^_.,„^(,_.) ^^^^^ ^^ _ ^3^^ 

dx J—oo Ztt 

For large yo (conjugate to x — *• 1), large perturbative corrections to the Fourier transform a{yQ, 6) 
exponentiate in much the same way as for the transverse momentum distribution in the Drell-Yan 
process discussed above, due to the factorization of soft gluons from the two energetic jets. As in 
the previous case, we neglect the nonleading contributions of coUinear gluons, and concentrate 
on the leading behavior in 6 only, given by 

cr(2/o, S) = exp{-S2j{yo, 6)) . (35) 

To lowest order the exponent S'2j is given by (compare (p!3D), 

S.,(y„, 5) = /C. / gj «,.(*) 2.SUk-) (-1^ - -1^) V - e""") . (36) 

with pi = ^(1, 1, 0) and p2 = ^(1, —1, 0) the light-like momenta of quarks in the two energetic 
jets. The function ^jet(^) is zero when fc is in a jet cone, and unity otherwise. Making explicit 
this kinematical restriction on the momenta of soft gluons, and including higher order corrections 
to S2j{yo,S), we get, to leading logarithmic accuracy, 

S2,iyo, S) = -2 /'' ^(l - e^'^o-) [^ ^T^^^AMk^)) , (37) 

where integration is performed over energy and over transverse momentum with respect to the 
axes of the outgoing jets. The soft gluon energy in (|37[) is restricted to be less than the energy 
of the "parent" quark, while the lower limit, kf > k^S, implies that the soft gluon is emmitted 
outside the jet cone. After substitution of (0) and ( p5| ) into (|3^ ) and integration over yQ, eq. ( pi] ) 
reproduces the leading perturbative series for the eikonalized 2-quark jet cross section with the 
total energy of gluons outside the cones given by (1 — x)Q. 

To take into account higher order corrections to the cross section, we have chosen, as usual, the 
argument of the coupling constant to be k^. An evolution equation for the exponent, 823 {yo, S), 
now follows from (pTj) (compare eq. ([I5|)) in terms of the "cone" parameter 6, 

^%1^ = 4 r ^(1 - e^'^y%^.,iaM6^)) . (38) 

dmo Jo ko ^ ' 

On the right-hand side of this equation, an infrared renormalon arises once again, from inte- 
grations over soft gluon energy ko ~ A/5, due to the singularity in ^^(/cq^^). Taking the one- 
loop expression (|^) for the cusp anomalous dimension and using the representation as(k\S^^ = 

10 



/q°° da{k'^5'^ / K^)""^^ , we may perform the integration over /cq in (|38|) to get an expression similar 
to (p!6D . However, an important difference from ([l6|) is that the IR renormalons now appear at 
a* = l/2/3i, l//?i, 3/2/3i, ... . Thus, in contrast with the transverse momentum distribution, the 
leading IR renormalon appears in the e"'"e~ —>■ 2 jets cross section at a* = l/2/3i, which gives rise 
to an 0{yoA/6) power correction to the exponent, S2j{yo, 5) in (|35|). We conclude that in the jet 
cross section, a{yo,S), the leading nonperturbative power corrections occur at level yoK/5. 

Now suppose we regulate the exponent S2i{iiQ,5) by some procedure, for example, a yo- 
dependent cutoff, by analogy to the Collins-Soper procedure in (|30|) and (0), or a principal 
value prescription for integrating over IR renormalons. We denote the result by 5'2j,pt(|/o; ^)- 
Taking into account the contribution of the leading IR renormalon to 5'2j(yo)'^)) we may write 
the cross section as 

a(|/o, S) = exp f-^2j,PT(l/o, S) - z^A2j + 0{yl)] , (39) 
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with A2J a new dimensionless nonpertubative parameter. In place of eq. (0), we now have 
d(T22j(x(5,(5) r^ dyi 



dx 



r- ^^_.,„Q(i_.+^,,A/g,) ^^p ^-S2,Myo, s)) . (40) 



As usual, A2J is dependent on our definition of S'2j,pt- We thus see that the effect of the leading 
infrared renormalon is to shift the scaling variable x by an amount inversely proportional to the 
typical transverse momentum in the jets, Q6. 

An operator intepretation of this result is not as simple as for the examples discussed in the 
previous section, because eq. (0) requires restrictions on both angles and energies of the soft 
gluons that give double logarithms. Nevertheless, the same leading logarithmic corrections can 
be generated by introducing a product of outgoing massive eikonal lines, defined as in eq. (p!0|), 
but with 

pl=pl = m\ m^ = Q^6^ . (41) 

The corresponding operator which takes into account soft gluon emissions into the final state is 
now 



U2j(x) 



$p,(oo,x)$], (cx),x) 



(42) 



and a corresponding cross section for the emission of gluons of total energy {1 — x)Q is given by 
(|3|) with 

a{yo,5) = We^e- = (0|4,t(0)f/2jct(2/)|0) , (43) 

where y'^ = (?/o,0,0). Indeed, eq. (^) has the same perturbative expansion as (^) and (p^), 
except that now the quarks have time-like momenta, pi = {E,p,0) and p2 = {,E,—p,0) with 
E^ = p^ + m'^, and the gluon momentum integrals are unrestricted. Then, integrating over k^ we 
get the following leading logarithmic expression for the exponent S2j{yo,S), 

52j(yo,5) = -2 rf^(l-e*'^o) /f ^§re..p(a.(A:n), (44) 

Jo Kq ^ ' Jklm?/E^ kf 

where E ^ Q/2 is the quark energy. Due to the cusp singularities of the Wilson lines, the 
/cq— integral in (|4^) is ultraviolet divergent and should be regularized, for example, dimension- 
ally, as in eq. ([T^) . The net effect of the regularization is to introduce a cut-off /x for large energies 
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ko. Taking fi = Q and identifying the mass of the quark as in (p|), we find that (|4^) coincides 
with expression (^) for S2]{yo,6). Clearly, to incorporate nonleading logarithms, which depend, 
in particular, on the internal structure of the jet, we must go beyond this rather crude approxi- 
mation. Besides these nonleading terms, there are additional corrections in this model associated 



with an ambiguity in the definition of the mass of eikonalized quarks. As discussed in [^ , the 
mass itself suffers from an ultraviolet renormalon, which introduces an ambiguity of order A into 
its definition. The corresponding ambiguity in the final-state energy contributes a term of order 
Ayo to the exponent exp[—iyoQ{l — x + A2]A/Q6)] in eq. (|40|) , down by a factor of 6 from the 
contribution we have identified. 

Since the first infrared renormalon singularity is proportional to yo, we may identify the 
parameter A2J with the first derivative of We+e-{yo,^) with respect to yo at yo = 0, by analogy 
to eq. (|25|) above, and we find that 



A 



2J 



(0|T 



%Moc){j^opM - 4pMHm] t[%MHM] |o)- 



(45) 



Note that this form involves only a single factor of the field strength, made gauge invariant by 
its combination with path-ordered exponentials. 

One may be tempted to suggest that such a matrix element, with only a single field-strength, 
must vanish. In fact, if we were allowed to remove the time and anti-time ordering from (^), 
this would be the case. That is, we could use the property $p$p = 1 of ordered exponentials 

to simplify (^5|) to (0| (j-'op^{0) — ^op^iO)) |0). This matrix element, however, is zero, because 
{0\^f,,p\0) = cp^, with c some constant, while p^(0|^^p|0) = (0|(p- D)$p<l>J|0) = from the 
definition of the Wilson line, which implies that c = 0. This argument fails, however, because 
the gauge fields on the pi and p2 eikonal lines do not commute in general. Thus, the matrix 
element in ( ^31) need not be zero. As usual, its precise value will depend on how we define the 
ambiguity in the resummed perturbation series. The important point here is that the matrix 
element depends on a single, although nonlocal, vacuum matrix element. A more detailed study 
of the large variety of jet cross sections will result in further information. 



5. Power corrections in inclusive lepton-pair production 

As final example, we consider the inclusive Drell-Yan cross section da-DviT, Q'^)/dQ'^, normalized 
to a structure function of deeply inelastic scattering, Fdis{x,Q^) [0, ^ ^, |2^, |13]. To be 
specific, the moments of the dilepton hard-scattering function are related to those of the DIS 
structure functions by 



Jo 



dri 



_id(TDY(r,Q^ 

dg2 



aoUj{n,Q'^) Fdis(?^,Q^) 



(46) 



where r = Q'^/s and o"o is the inclusive Born cross section, as in eq. (|1]), and where -F'dis(^)Q^) 
is the corresponding moment of F-£,is{x , Q"^) with respect to x. For large values of n (conjugate 
to r — > 1 and a; — *■ 1), iu{n, Q"^) contains large perturbative corrections of the form aj In n. 

The logarithms of n, however, are readily resummed. In fact, up to corrections that are 
bounded in n, the moments of the hard-scattering function exponentiate in the form 



Lo{n,Q^) = RiasiQ^)) exp E{n,Q 



(47) 
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where R{as) is a finite function of as, wliile tlie function E{n,Q'^) generates all logarithmic 
corrections in n from the integrals^ 



E{n,Q^ 



/■I , z 

-2 / dz — 

io 1 



n— 1 



-1 



L r ( 

2 -■- cusp I 



(1-^)2Q2 /c^^ 



a,(A;2)) + r(a,((l-^)Q' 



(48) 



The anomalous dimensions Fcusp and F have been encountered above in eq. (P). E{n, Q"^) contains 
infrared renormalon singularities from small values of kf. We may isolate the leading infrared 
renormalon by interchanging the kf and z integrals in the first term of (^8]), 



E{n,Q' 



^ I t-2 cusp 



/O n,t 

2(n- 1 



Jl~-kt/Q 










cusp 



«.(A:D) 



I- z 
kx k"^ fn — 2 



4 



+ 



(49) 



where we have expanded z"~^ in powers of 1 — z, and have suppressed (1 — z)^ and higher. 
Such terms produce higher powers of nkt/Q. For large values of kf, all such terms contribute to 
the perturbative expansion of E{n,Q'^). The leading infrared renormalon, however, which may 
be identified in the same manner as above in eq. (|TBp, is present only in the first term in the 



expansion of (^91) , which we denote 

4nADY(g) 

Q 



n 



= 4-- / d/cT Fcusp(as(fcT)) • 
Q Jo 



(50) 



Again, Ady{Q) includes both perturbative and nonperturbative contributions. The latter is 
precisely the 1/Q contribution identified in [T^|. Following our procedure above, we now find a 
nonlocal matrix element which corresponds to Ady{Q)- 

We begin by defining a product of Wilson lines that reproduces the logarithmic behavior of 
the Drell-Yan cross section near the boundary of phase space ||24|| , 

dyo 



WBYiko,fi) 



27r 



-^'°y^{o\uUo)uMym 



(51) 



with Udy the product of ordered exponentials in the directions of the incoming particles, given 
above in eq. (|10|) ioi x = y = {yo,0). As indicated, VFdy requires renormalization at the scale 
fj,. From the procedure introduced in [^], Wdy satisfies an evolution equation whose solution is 
the convolution form^. 



WBY{ko,fi) 



oo Qn n „i 

R{as,e)Y^ —rY[dzi 



„=o ^! ^Vl -^0 



Zi Jo 



(i-Zi)/i dk 



t,i 



K 



Tr 



2+2e ^ cusp 



{O^sikl)) 



+ t{as{{l~z)ii^ 



5 [h-Y^i.^^- Zi)^i 



(52) 



Here, R is an (infrared sensitive) function of a^ and e, and F organizes nonleading logarithms, 
which may differ from those in E{n, Q'^) above. We may now evaluate the Fourier transform of 
Wdy explicitly. 



Wdy(i/o,/^) 



dk 



^ Jkoyo 



27r 



e^'=°^"iyDY(fco,/x) 



(53) 



^In |£^ and |g^, 2rcusp and 2T are denoted gi and 52, respectively. 
^The reasoning is essentially identical to that given in Section 5 of |21[ 
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We then recognize that its derivative with respect to yo, evaluated at yo = 0, generates the 
product of Wdy{0, fi), which can be normahzed to unity, and a single integral over kt, in which 
appears in the scale of the running coupling, 



.dlnWoYiyo, fj.) . ^ f^"' dk 



2 



-^ "^^"""1 ,0=0 = 2/- -i^r>,(A;,2))(/i-fc,). (54) 

oyo JO ft J 

While both of the terms on the right-hand side contain an infrared renormalon, only the second 
corresponds to the function ADY{fi) identified from the resummed cross section above. To isolate 
^dy(/^) from the ultraviolet contribution in eq. (0), we simply form the combination, 

. ._ i dlnWBY{yo,fi), 1 92lniyDY(z/o,/i)| .^.^.x 

Ady(/x) - ^^^ 1,0=0 - - ^^2 1.0=0- (55) 



Here, the second term, which is evaluated just as in (|5^), serves only to cancel the large, ultra- 
violet contribution. Its nonperturbative contribution, however, begins only at order /i~^ ~ Q~^, 
which we neglect in our approximation. Note that the second term is not unique, and higher 
derivatives (times higher powers of 1/yu) would equally well cancel the ultraviolet contribution, 
while contributing infrared renormalons that are yet further suppressed in /x ~ Q. 

6. Interpretation 

We have studied several cross sections, whose leading logarithms are generated by products of 
ordered exponentials. In each case, we found that nonleading power corrections are required 
by the presence of infrared renormalons in the corresponding resummation formulas. These 
corrections, in turn, may be represented as vacuum matrix elements of field strengths, integrated 
over the paths of the original ordered exponentials. It is important to stress that our results 
do not depend on an extrapolation of the k^ integral in, for example, eq. (^), to soft momenta, 
k^ ~ A^. Rather, eq. (^ organizes the large-order behavior of the cross section that follows from 
the evolution equation (|T3]). In this sense, the presence of the divergence in the running coupling 
follows from the factorial behavior of perturbation theory, not the other way around. 

In the specific cross sections we have studied, we have found that the first nonleading-Q 
behavior is associated with nonperturbative matrix elements, such as (pSf), ( ^5]) and (0). From 
our critereon of consistency for the full theory, including both perturbative and nonperturbative 
contributions, we conclude that these power corrections exponentiate, to give, for example, a 
nonperturbative Gaussian distribution in b. The Fourier transform of such a function gives a 
Gaussian behavior in momentum space, which must be convoluted with the perturbative Qt 
distribution. On the other hand, behavior linear in the transform variable, such as n/Q in 
the Drell-Yan normalization and yo/Q in the jet cross section, is associated with a shift in the 
conjugate kinematic variable. 

The reasoning that leads to this picture is similar to that described in [Q, which connected 
infrared renormalons in the total e"'"e~ annihilation cross section with the local condensate 
(0|F^(0)|0). Here, we must go beyond the class of local operators to "nonlocal condensates", 
but with the benefit of greater flexibility, and perhaps, generality of application. Of course, this 
procedure lacks the guiding principles of the operator product expansion, and we have chosen 
to let perturbation theory suggest for us the form of nonperturbative structures that we may 
expect. While other contributions, not related to perturbation theory, are probably present, 
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the self-consistency of the theory demands the presence of those we identify from perturbative 
calculations. 

We emphasize again the necessarily ambiguous nature of the magnitudes of the higher-twist 
matrix elements (e.g., (|25|)), which depend on the manner in which the perturbative integrals 
have been constructed and, indeed, the order to which they have been computed. In some 
sense, however, the situation is a bit better than for the total e"'"e~ annihlation cross section, as 
described in Ref. ||2^, where beyond low orders, the perturbative series will begin to show T{n) 



behavior. In resummation formulas, such behavior is already included with known coefficients. If 
we succeed, therefore, in constructing a perturbative exponent S whose (arbitrary) higher twist 
is not overly sensitive to higher-order corrections to (for instance) Fcusp, then the value of our 
nonperturbative matrix elements will correspondingly be stable to higher-order corrections in the 
perturbative calculation. Of course, there are other short-distance corrections not included in 
the resummation, but these will appear only at still higher twist. 

Our discussion has been quite formal, but we hope that it will prove valuable to recognize 
that a limited set of nonperturbative parameters may enter into the first higher-twist corrections 
to many high energy cross sections. Such "universality" has proved of great practical use when 
applied to local condensates in QCD sum rules, despite the limitations and ambiguities inherent 
in their combination with perturbative calculations. Beyond this, were it possible to actually 
compute, by some nonperturbative method, highly nonlocal and relativistic operator combina- 
tions such as eq. (^), in a manner consistent with a particular perturbative construction of S, 
then the formalism presented here would provide new tests of the theory. 
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